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ABSTRACT 

The CQUAD4 thin plate element implemented in COSMIC NASTRAN is 
capable o5 Q modeling thin layered plate and shell structures composed 
of orthotropic lamina. Fiber-reinforced composites are *** 

classes of inhomogeneous and non-isotropic materials which can be 
treated Although the CQUAD4 has been extensively checked in static 
cases little validation has been carried out for vibration response 
modeling This paper documents validation of the CQUAD4 element s 
accuracy for vibration response analysis of thin laminated composite 

plates. 

The lower-order natural freguencies and mode shapes of ten 
glass fiber-reinforced plastic (GFRP) and carbon fiber-reinforced 
clastic (CFRP) plates are computed and compared to published expen 
menta^and numerically-computed data. A range of ply geometries 
including unidirectional, cross-ply, and angle-ply are considered. 
The plates' length-to-thickness ratios all lie in the vicinity of 
100 to 150. The CQUAD4 plate idealizations provide natural frequen- 
cy predictions within ten percent of measured data for all six 
lowlst modes of seven of ten pistes. For tw ° ? £ . temsining three 

plates, only the fundamental frequency is predicted with an er 
qreater than ten percent. Results for the one remaining plate do 
not correlate with published data, possibly because e *roneous 
reporting of its geometry or material properties in the literature. 
To P obtain accurate frequency predictions, lamina in-plane elastic 
moduli had to be tuned to reflect each plate's fiber volume 
fraction. 

These results show that the NASTRAN CQUAD4 plate element is 
useful and reasonably accurate for vibration and shock analysis of 
structures composed of thin fiber-reinforced plastic plates. 

INTRODUCTION 

There is strong current Navy interest in exploitation of fiber- 
reinforced plastics as lightweight materials for a wide variety of 
ship structures. These structures must be designed to withstand in 
service loads of quasi-static, transient dynamic, and steady-state 
dvnamic nature. For many ship structures, transient shock is a pn- 
S ioaS? Fo; these and others, steady-state vibration response 
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t J} e 1 shi P ,s acoustic signature. In both cases, the struc- 

are 

e employed to obtain realistic modal property estimations. 

The NASTRAN finite element code is one of the Navy's premier 
toois for steady-state vibration response analysis of ship and 
submarine structures. Undamped natural mode analysis and forced 

bi NAsi^ N re aI°i;eli a 2s 1 nodf] W f th hysteretic damping can be performed 
£ NAbiKAN, as well as modal frequency and loss factor analvsis for 

structures with viscoelastic damping materials (Ref. 1 ) . NASTRAN is 

?ad?J.S ey C0, " p0n ? nt ° f the NASHUA suite codes for performl^ 

ged JtouctotoS a Ref aC °“^° scatteri "9 analysis of vibrating submer- 

. . COSMIC NASTRAN CQUAD4 element is designed to model aniso- 

The P th. layere a plates as well as homogeneous and isotropic plates 
The theory and assumptions behind the CQUAD4 have been informally 
documented (Ref 4) . The CQUAD4 has been found to be more Curate 
for a given finite element grid than its predecessor, toe cSSam 
element, for prediction of low frequency eigenmodes of thin-walled 
cylindrical shells composed of isotropic materials fRef Tn 

* n0 r^ d9 % n0 comparable study of the ^u^y'of’the 
CQUAD4 formulation for anisotropic plate or shell vibration vet 
exists, particularly for plates and shells composed of of 

fiber-reinforced plastic lamina. yers ot 

This paper summarizes an investigation of the accurarv n f 
NASTRAN CQUAD4 membrane and plate bending element for vibration 
analysis of structures composed of thin fiber-reinforced composite 
plates. This is accomplished by comparisons of NASTRAN-computed 
undamped natural frequencies of ten GFRP and CFRP plates Sit h 
published experimental data and other numerical predictions. 

NOTATION 


E 11 

Lamina extensional modulus in 
direction parallel with fibers 

E 22 

Lamina extensional modulus in 
direction transverse to fibers 

E 33 

Through-thickness extensional 
modulus 

G 12 

In-plane lamina shear modulus 

G 13 

Transverse lamina shear modulus for 
out-of-plane shearing of a fiber 
cross section 


52 



G 2 3 

Transverse lamina shear modulus for 
in-plane shearing of a fiber cross 
section 

L 

Plate side length 

t 

Plate thickness 

V f 

Fiber volume fraction 

V 12 

In-plane shear strain 

v V 

• 13' 23 

Out— of— plane shear strains 

£ 1' € 2 

In-plane extensional strains 

e 

Fiber orientation angle with 
respect to one plate side 

^ 12 ' ^21 

In-plane lamina Poisson's ratios 

^ 23 

Transverse lamina Poisson's ratio 

<? 

Average mass density 

^1' ^2 

In-plane extensional stresses 

1 12 

In-plane shear stress 

'Y *Y 

' 13' 23 

Out-of-plane shear stresses 


ABBREVIATIONS 

avg. 

Average 

CFRP 

Carbon fiber-reinforced plastic 

CLT 

Classical lamination theory 

consis. 

Consistent 

CPT 

Classical plate theory 

DOF 

Degree (s) -of-freedom 

FEA 

Finite Element Analysis 

GPa 

Gigapascals 

GFRP 

Glass fiber-reinforced plastic 

Hz 

Hertz 

kg 

Kilograms 
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m 

Meters 

mm 

Millimeters 

params . 

Parameters 

psi 

Pounds (force) per square inch 

ref. 

Reference (s) 

RMS 

Root-mean-square 

sym 

Midplane-symmetric 


dACJmjKOUND 


Anisotropic plate analysis is more difficult and involves more 
variables and parameters than isotropic plate analysis. The impor- 
tance of accounting for transverse shear flexibility in relatively 
thin composite laminates is discussed in the next section. The ways 

described ° 4 addresses these difficulties are also briefly 


TRANSVERSE SHEAR DEFORMATION ISSUES 

In composite plate mechanics, the counterpart of the wel l -lfnnun 
Kirchhoff (Classical plate) theor^ (CPT) for isoSopL plates is ?he 
Classical Lamination Theory" (CLT) . The two theories 

iEI J he ^ e klnematlc assumptions regarding the deformation of 
the plate with respect to its middle surface; that is, sections 
originally planar and perpendicular to the middle surface remain 
planar and perpendicular in the deformed state. The mathematical 
evelopment of CLT occupies much of the text by Jones (Ref. 6). The 
reader should consult this (or some other) text for detailed 
exposition of CLT assumptions and derivation of the CLT equations. 

nla . The accuracy of the Kirchhoff kinematic assumption in isotropic 
plate mechanics degrades when the plate thickness becomes signifi- 
cant compared to its span length. For isotropic metallic plates 
transverse shear-stiffness-governed transverse deflection becomes 
significant: f'? latlve to flexural deflection when the span length-to- 

^?^ neSS ratl ° (L/ ^ is suffi ciently small. An idea of required 
smallness can be gained from the discussion of "corrected" plate 
flexural waves in Chapter II, Section 3b of Cremer, Heckl and 

an 9 ?L^ 6f: 7) \ The ^ Sh ° W ' f ° r a uniforinl y thick plate composed of 
an isotropic material, that transverse shear effects decrease flex- 

ural wavespeecl by ten percent when the flexural wavelength is about 
six times the plate thickness. For one-half wavelength over a plate 
span length, this limitation translates to an L/t ratio of 3. 

1 , , ' n ,^ e of high modulus composite plates, an analogous 

limitation of the adequacy of CLT is encountered at large r L/t 
This more stringent limit arises because the ratio of effective 
laminate extensional elastic modulus to shear modulus is a key fac- 
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tor (other than section geometry) governing the magnitude of shear 
deformation relative to flexure. A rough idea of limiting L/t 
ratios for the adequacy of CLT for composite plates follows from 
insertion of some "ballpark" ratios of effective laminate extension- 
al and transverse shear moduli for GFRP and CFRP into the approxi- 

^expression for shear-corrected plate flexural wavespeed found 

on page 115 of Reference 7. Using Ej- / G - i3 for GFRP and 
forCFRP ten percent differences mHavespsed arise for L/t of 6.75 
and 20 . 9 ' respectively. (A comparable modulus ratio for isotropic 
materials is 2 6) . These ad hoc assessments are qualitatively corro 
SorSd by sta^c examples found in section 6.5 of Reference 6 In_ 
the problem of cylindrical bending of a CFRP strip with E / G 
12 I, maximum static deflection predicted by CLT 1 . twenty jbercefrt 
smaller than the true shear-corrected solution at L/t - 20. 

The low frequency composite plate vibration literature is domi- 
nated by evaluation of methods for account of transverse (^terlami 
nar) shear in prediction of natural vibration frequencies and modal 
deflections. The inadequacies of CLT even for fundamental plate 
frequencies are repeatedly demonstrated in the literature for L/t 
frequency Approaches have been developed to provide finite 

element-based plate formulations to handle through-thickness stress 
fields for arbitrary L/t. In these complex formulations, transverse 
normal stresses are no longer assumed to be zero, transverse shear 
stresses are constrained to be zero on upper and lower laminate 
surfaces, and shear stress continuity between laminae is maintained. 
As a result, sections perpendicular to the plate midplane rotate 
with respect to the midplane and warp out of a planar configuration. 

The plate finite elements simulating through-thickness stress 
fields in laminates are, obviously, quite sophisticated. Tables II 
and VI of Mallikar juna and Kant (Ref. 8) provide a good flavor for 
the performance of some higher-order approaches for reckoning with 
the effect of transverse shear and normal strains in plate 
freouencv prediction. Table III shows that CLT provides reasonable 
fundamental frequency predictions for a simply-supported CFRP angle 
ply^plate for L/t at and above 20. Even for such highly anisotropic 
plates, the complications of very complex high-order transverse 
shear theories, necessary for laminate strength and structura 
integrity°probl ems , are not justified in vibration problems unless 
L/t is less than about 20, or if short-wavelength (high frequency) 
vibration modes are of interest. 

The NASTRAN CQUAD4 membrane and bending element embodies the 
assumptions of CLT, but contains first-order corrections for trans 
verse P shear flexibility. When the CQUAD4 is used to 
eous plates, transverse shear strains vary linearly with the thick 
ness coordinate, and zero shear stress and strain boundary condi- 
tions on upper and lower plate surfaces are not satisfied. The 
shear energy implied by the linear distribution is corrected by a 
multiplicative constant to produce the energy implied by the true 
quadratic distribution. As a result of these assumptions, sections 
perpendicular to the plate midplane are allowed to rotate with 
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respect to the midplane, but remain straight when the plate is in 
the deflected state. This ad hoc transverse shear correction ex- 
tends the element's range of validity to thicker homogeneous plates 
but is not acceptable for inhomogeneous layered plates. 

When modeling inhomogeneous plates composed of orthotropic 
layers, a quadratic transverse shear strain distribution is assumed 
in each layer of the CQUAD4 element. Interlaminar shear strains are 
matched at lamina interfaces, and zero shear stress boundary condi- 
tions are enforced on the upper- and lower-most lamina surfaces. 
However, the through-thickness normal stress is assumed to be zero 
and complete consistency between the strain-displacement equations' 
for m-plane direct strain and transverse shear strain is not main- 
tamed. ^ The CQUAD4 is thus seen to overcome the limitations of CLT 
for laminated plates, but does not represent all aspects of the 
kinematics of three-dimensional continua taking place in relatively 
thick laminates. 1 

The CQUAD4 element is discussed in more depth in the following 
section. ^ 


THE CQUAD4 ELEMENT 


The CQUAD4 is a four-noded planar element possessing membrane 
flexural, and transverse shear stiffness. In the case of layered 
plates, individual laminae are not modeled explicitly; rather 
equivalent stiffness matrices for the plate as a whole are defined 
Each lamina is assumed to be in a state of plane stress, and the 
laminae are presumed to be perfectly bonded by infinitesimally thin 
non-shear-deformable layers. Each lamina is assumed to be specially 
orthotropic, with six independent elastic moduli when through- 
thickness direct stresses and strains are ignored. Any alignment of 
lamina fiber axis with respect to the local element coordinate 
system can be accomodated. Hence, any layup or stacking sequence 
can be handled (unidirectional, cross-ply, regular or irregular 
angle-ply) . Layups unsymmetric with respect to plate midplane can 
also be modeled, as membrane-bending coupling is accounted for when 
it occurs. Each lamina may also be composed of a different ortho- 
tropic material, if desired. 


The element's stiffness matrix terms for determining in-plane 
displacements and flexural rotations as a function of imposed forces 
and moments arise from CLT assumptions. The CQUAD4's force-versus 
f^ ra11 } equations for membrane, flexural, and membrane-flexure coup- 
ling (Ref. 4) are identical to those developed in sections 2 1 
through 2.6 of Reference 6. The kinematic assumptions regarding in- 
plane and flexural strains and displacements follow classical 
assumptions and require no explanatory remarks here. Reference 4 
documents the force-versus strain matrix terms for transverse shear 
strains. These are based on overall element equilibrium, continuity 
of transverse shear between adjacent laminae, and satisfaction of 
zero shear strain and stress boundary conditions on the upper- and 
lower-most laminate facings. As mentioned earlier, the strain- 
versus displacement matrix is based on the assumption that through- 
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thickness planar sections rotate with respect to the plate raidsur 
face and also distort out of planes originally perpendicular to the 

midsurf ace . 

The CQUAD4 is an isoparametric element, whose displacement 
fields are interpolated through space by linear variations of 
plane and transverse displacements (and rotations about midsurface) 
between grid points. The associated in-plane strains are constant 
between grid points but vary linearly with the thickness coordinate. 
However, transverse shear strain varies quadratically through each 
lamina thickness. 

Those who use the CQUAD4 to model fiber-reinforced composite 
plates must realize that the input elastic moduli are effectiv e 
moduli for a particular fiber -matrix combination- There are many 
different fibers in use (glass, carbon, kevlar, and boron are exam- 
ples) and many resins or matrix materials, each of which has their^ 
own unique elastic moduli. Although matrix resins are usually con 
sidered to be isotropic, fibers have distinct extensional, trans 
verse, and shear moduli. The effective in-plane moduli of an ortho 
tropic continuum defined to be equivalent to the actual inhomogen- 
eous fiber and resin system are (sometimes nonlinear) functions of 
the extensional, transverse, and shear moduli of the fiber, the 
extensional and shear moduli of the resin, and the Poisson's ratios 
of the fiber and resin. The fraction of the lamina volume occupied 
bv fiber material is an important variable defining the magnitude of 
effictiSe in-plane lamina moduli. Some strength-of-materials rela- 
tionships defining effective lamina in-plane moduli as a function of 
constituent matrix and fiber moduli and fiber volume fraction are 
developed in sections 3.1 and 3.2 of reference 6. Reference 9 pro 
vides a handy tabulation of effective lamina elastic moduli under 
the assumption of a transversely isotropic lamina. Chamis provides 
formulas for effective out-of-plane shear moduli as well as the more 
commonly reported in-plane extension and shear moduli. 

APPROACH 

Specifics of the present vibration modeling study are now 
described. 

PLATE VIBRATION SPECIMENS 

It was desired that the validity of the CQUAD4 be proven by 
comparison to numerical results obtained independently by other 
researchers, and to experimental data, if possible. Lin, Ni, and 
Adams (Ref. 10,11) and Xiao, Lin, and Ju (Ref. 12) have published 
experimental data and finite element computations for the lowest six 
vibration modes of square CFRP and GFRP plates for fre «' . 

strained boundary conditions. They measured and computed both natu 
ral frequencies and damping loss factors, which makes their work 
almost uniquely complete and thorough. Reference 11 contains data 
for nine plates repeating that for four plates treated in reference 
10. Reference 12 revisits three previously examined plates with 
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more refined finite element formulation allowing plane sections to 
rotate and warp, with consistent correction of all strain-displace- 
ment relations for these kinematic conditions. This paper contains 
data for one additional new plate, making a total of ten unique 
plates in all (four CFRP and six GFRP) . Refinement of their finite 
element method (Ref. 12) improved correlation of their predictions 
with measurements. 


The geometric and material parameters of the ten plates studied 
by Lm et. al. and Xiao et. al. are listed (in SI units) in Table 1. 
The parameters given in reference 12 for plates 770 and 772 are 
inconsistent with those reported in reference 11. The NASTRAN study 
confirmed that the side lengths, thicknesses, average mass 
densities, and fiber volume fractions given for plates 770 and 772 
in reference 11 are the correct values. Further correlation of 
Pl a ^ e parameters, frequencies, and mode shapes between references 
10-12 revealed accidental reversal of mode shape plot labels in ref- 
erence 11. in addition, two unidirectional GFRP plates of different 
size are reported in references 11 and 12 with the identification 
number 761. (They are herein distinguished from the other as 761L 
and 761X) . These discrepancies initially caused much confusion, but 
the author is confident that the data in Table 1 is correct. 


Some features of these plate specimens are notable. All of 
them are square, and have L/t ratios in the vicinity of 100 to 150. 
Even though such L/t would seem to be in the range of applicability 
of CLT, Xiao et. al. show that CLT, which totally ignores transverse 
shear deformations, overestimates natural frequencies by factors as 
high as sixteen percent over a theory with first-order shear correc- 
tion. All of the laminates listed in Table 1 are symmetric about 
the plate midplane, eliminating flexure/extension coupling effects 

1 ™* ( ^ FRP s P ecimens 734 , 761 ( L) , 761 (X) and CFRP specimens* 

762 and 764) have "specially orthotropic" lamina (all fibers aliqned 
with the plate sides). For these plates, there is no coupling 
between in-plane extension and in-plane shearing, so their vibration 
mode shapes have nodal lines more or less parallel with the plate 
edges. The other five plates (765, 769, 771, 770, 772) have at 
least some plies with fibers angled relative to plate edges pro- 
viding more complex mode shapes. ' 


All plates were tested with "free" edges (supported by soft 
foam rubber strips) and numerically analyzed with zero-constraint 
boundary conditions. (Although many practical design problems of 
interest would involve plate structures with boundary constraints 
this feature eliminates uncertainties about which boundary degrees- 
of-freedom (DOF) to constrain to obtain "simply supported" edges) 

An iterative technique was used to obtain natural frequencies and 
modal damping loss factors. First, the specimens were excited into 
steady-state vibration by an electrodynamic shaker, and approximate 
resonance frequencies were determined. Nodal lines for each excited 
natural mode were located by the classical Chladni sand pattern 
technique. Then, for each mode of interest, locations of the 
support strips were then adjusted to align with nodal lines. A 
transient excitation technique was then used to obtain more precise 
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estimates of resonance frequencies and modal damping loss factors. 
LAMINA MATERIAL PROPERTIES 

Effective in-plane lamina elastic moduli are reported in Lin 
et. al. and Xiao et. al. These are listed in Table 2, in both SI 
and English units. These moduli pertain to microstructure 
where half of the total ply volume is occupied by fibers and half by 
the matrix. For each plate analyzed here, the moduli input to 
NASTRAN must be adjusted up or down from these nominal values, 
according to the measured fiber volume fraction (V ) for each plate : 
In general, the moduli vary nonl inear ly with V f , bfit experimentally 
verified semi-empirical equations are available for determining V~ 
adjusted moduli. References 10-12 reported only the nominal moduli 
in Table 2, but provided a literature source (Ref. 13) for accom 
plishing the adjustments. 

For orthotropic laminae in a state of plane stress (in the 1-2 
plane), the only independent engineering moduli are E-w E, 2 , G 12 , 
and 'v . The additional Poisson's ratio V 21 must satisfy the 

relationship: 


12 


/ E 


11 


^21 / E 22 


The NASTRAN CQUAD4 element enforces this constraint on (Ref* 

4) . With these moduli, all in-plane strain and stress componen s 
are defined by the CLT. However, the NASTRAN CQUAD4 element also 
requires, as input, nonzero transverse shear moduli G _ and G , 
associated with out-of-plane shearing, which are not Specified or 
required in CLT. Fortunately, fiber-reinforced lamina may often be 
assumed to be "transversely isotropic" in analysis of composite 
structures, as implied by reference 9. This means that if the 
lamina lies in the 1-2 plane, and the fibers are aligned in the 1 
direction, then the transverse shear modulus G-- is equal to the 
in-plane shear modulus G. _ . The transverse shSar modulus associated 
with shearing of the matrix "around" the fibers and distortion of 
the fiber in a plane perpendicular to its axis (G 23 ) remains to be 
determined. To clarify these matters, the stress and strain compo- 
nents for a transversely isotropic lamina with zero direct stress 
normal to the lamina are illustrated in Figure 1. 

Lin et. al. and Xiao et. al. did not report either of the 
transverse shear moduli used in their analyses. Educated guesses 
had to be made for both the CFRP and GFRP G . Reference 8 provides 
nondimens ional elastic moduli for a CFRP-lixe material which J*?, 
closely satisfied by the in-plane parameters for the CFRP in Table 
2. According to reference 8, G_ 3 should equal 0.2E 2 _. This implie 
value of G„ was used in the prSient CFRP plate analyses, based on 
fiber volu&l fraction-adjusted E For the glass fiber lamina, G 23 

= 0 . 6 G was assumed, based on in-house experience with such mater- 

ials. 1 ^he assumed out-of-plane moduli are indicated m the notes 
section of Table 2. G 23 could also be estimated by methods reported 

in reference 9. 
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Table 3 contains lamina elastic moduli for each of the ten 
plates, after adjustment of the nominal moduli in Table 2 for fiber 
volume fraction. Lin et. al. did not correct the in-plane Poisson's 
ratio for fiber volume fraction, even though Ni and Adams show that 
it decreases appreciably for increasing V . Lin et. al. claim that 
accounting for this decrease had no significant effect on their 
calculated natural vibration frequencies. We match their assumption 
by keeping constant at 0.3 for all ten plates, in general 

V 12 should be adjusted for V f . 

ASSUMPTIONS REGARDING LAMINATE BEHAVIOR 

The kinematic assumptions of the CQUAD4 element have been dis- 
cussed previously; namely, that planar sections perpendicular to the 
plate midplane can rotate and warp relative to the midplane. NAS- 
TRAN-computed natural frequencies for CQUAD4 plate idealizations 
will, subsequently, be compared to some CLT predictions (Ref. 12) 
and to the FEA predictions of Lin et. al. (Ref. 11) and Xiaoet. al 
(Ref. 12) These three approaches differ in accuracy, and it is 
important to understand how the NASTRAN predictions should compare 
with them. Kinematic assumptions are discussed first. 

As discussed previously, CLT totally ignores any transverse 
shear deformation effect, and will always predict lower— order 
^^■krat ion frequencies which are too high for laminates with low L/t 
and high extensional-to-shear modulus ratios. In terms of through- 
thickness lamina kinematics, CLT prescribes linear variations of 
direct and in-plane shear strains, and zero transverse shear 
strains. The kinematic assumptions of Lin et. al. (Ref. 11) and 

C QV AD4 ( Ref * 4) are identical, and imply linear through- 
thickness variations of direct and in-plane shear strains and 
quadratic through— thickness variations of transverse shear strains 
with zero shear stress boundary conditions on the upper— and lower- 
most lamina facings satisfied. However, this thick plate-type 
theory is only an approximation; consistent correction of all 
strain-displacement relationships when cross-section warping is 
allowed requires cubic through-thickness variations of direct and 
in-plane shear strains. Xiao et. al. (Ref. 12) include this consi- 
derable complication in their plate-type FEA formulation, which can 
be understood as a special case of three-dimensional elasticity. 

From an understanding of kinematic assumptions alone, CQUAD4 
composite plate idealizations should provide natural frequency 
estimates that are (1), more accurate than CLT, (2) as accurate as 
the Lin et. al. predictions, and (3) less accurate than the Xiao et 
al. predictions. However, other factors will influence the com- 
parisons between FEA-predicted frequencies; particularly, the 
polynomial form of the interpolation functions expressing element 
displacement fields in terms of grid point displacements, and the 
mass matrix formulation (consistent or lumped) . The CQUAD4 utilizes 
linear interpolation functions with respect to the element's four 
grid points. Lin et. al. do not specify the interpolation func- 
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tions used in their "8-node 4 O-degree-of -freedom" elements, although 
they can be guessed to be of at least quadratic order. 

NASTRAN ANALYSES 

The results of the NASTRAN computations are now compared to 
the calculations and experimental data in references 10-12 . 

FEATURES OF THE NASTRAN PLATE IDEALIZATIONS 

Lin et. al. and Xiao et. al. employed a 6-by-6 mesh of 8-noded 
isoparametric rectangular elements in all of their plate FEA ideals 
zations. Each element possessed 40 DOF, (5 per grid point), imply 
inq that all rotations about axes normal to the plate surface were 
constrained. As in the CQUAD4 , they applied a numerical condition- 
ing factor to transverse shear stiffness terms to eliminate excess 
ive shear stiffness ("shear locking"); a consequence of numerical 
integration of element stiffness. The NASTRAN CQUAD4 elements also 
have five DOF per node, with normal axis rotation constraint pp 
and are also conditioned to avoid shear locking. 

The CQUAD4 differs from the Lin and Xiao et. al. elements m 
one important way; they are 4-noded isoparametric quadrilaterals and 
thus have linear interpolation of in-plane displacements between 
arid pain+’ c instead of quadratic interpolation. Thus, although t 
Si et al element is kinematically similar to the CQUAD4 as far as 
through-thickness shear effects are concerned the 

in-plane displacement fields differ as a function of the plate s in 
plane dimensions. The Xiao et. al. element provides for a more com- 
plex through-thickness displacement and strain distribution than the 
Lin et. al. element and the CQUAD4 , as discussed earlier. 

The NASTRAN CQUAD4 meshes used here consist of a 12-by-12 grid 
of elements, with 169 grid points and 845 unconstrained DOF. These 
idealizations are roughly comparable to the Lin and Xiao models in 
modal displacement field interpolation quality. Three of the plates 
were initially modeled with 6-by-6 meshes with 49 grid points and 
245 DOF, but these idealizations did not provide sufficient mode 
shape resolution to be acceptable. The Lin and Xiao et. al. mesh is 
compared to the NASTRAN mesh in Figure 2. 

The edqes of the NASTRAN plate idealizations were unconstrain- 
ed and all rotations about axes normal to the plate surface were 
suppressed. A SUPORT input record imposing fictitious constraints 
on all five remaining DOF on one grid point was used to provide 
zero- frequency rigid body modes. The FEER eigenmode extraction 
method was used in NASTRAN Rigid Format 3 (normal modes analysis) , 
with the lowest twenty modes requested. The Inverse Power method 
was employed in some trial runs; it provided about half as many 
frequencies as FEER but at more than two times greater run time and 
cost. All computations were performed by RPK COSMIC NASTRAN, 1990 
release, installed on the DTRC CRAY X-MP supercomputer in a COS 
operating system environment. Typically, 22 to 24 modes were ex 
tracted by FEER in 23 to 24 CP seconds. 
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Lin and Xiao do not mention whether they used lumped or con- 
sistent mass matrix formulations. Both options were considered in 
the NASTRAN study. 

NASTRAN MODAL ANALYSIS RESULTS 


NASTRAN analyses utilized nominal elastic moduli based 
on fifty percent fiber volume fraction. Although predicted frequen- 
cies were in the general vicinity of measured values, they did not 
compare with CLT results, the Lin et. al. results, or the Xiao et. 
al. results in the expected way. That is, the CQUAD4 did not neces- 
sarily appear more accurate than CLT and less accurate than the Xiao 
et. al. approach. It suffices to say that fiber volume fraction 
strongly influences effective moduli and must be accounted for to 
obtain credible composite plate modal predictions. 


Relative NASTRAN frequency predictions are compared with CLT 
Lin's, and Xiao's results in Table 4. There is, in general, good' 
correlation with measurements, and NASTRAN 's frequency error trends 
are comparable to those in the Lin et. al. analyses (especially 
those for a consistent mass matrix). The fact that effective lamina 
elastic moduli must be corrected for fiber volume fraction to obtain 
credible natural frequency predictions for laminated composite 
« a ^ GS 1S em P* las i 2 ed in Table 5, where root-mean-square (RMS) values 
of frequency prediction percentage error are tabulated for the six 
modes of each plate. RMS error is seen to be significantly reduced 
m most cases when fiber volume fraction-corrected lamina moduli are 
employed. The major exception is plate 765, which still suffers 
from some large unknown systematic error making predicted frequen- 
cies far too low. Results for plate 770 were not much changed since 
actual V was already close to one-half. RMS errors for plates 764 
and 771 remain at ten percent and above. However, the major part of 
these high errors involves their fundamental modes, which can be 
easily impacted by test article boundary constraint. (The Lin et. 
al. analysis also predicted overly high frequencies for these 
modes) . 


• ^ ^ bS m 1 ^! :e measu ^ ed a ? d predicted natural frequencies are summar- 

ized in Table 6, which lists Lin et. al. measurements and computa- 
tions and the NASTRAN CQUAD4 idealization results for lumped and 
consistent mass. Mode-by-mode and average percentage frequency 
prediction errors are listed in Table 7 along with mode shape de- 
scriptions. The NASTRAN lumped mass model gives the lowest averaqe 
error in seven of ten cases, and the NASTRAN consistent mass model 
is best m two of the three remaining cases. (Interestingly the 
consistent mass model errors roughly parallel those of the Lin 
model). Lin et. al. calculations are ''best" only for plate 765 for 
which the author believes a parameter was misdocumented in reference 
11. In six of the nine plates other than 765, NASTRAN predicts 
frequencies that are all at or within ten percent of measured. In 
the other three, errors are greater than ten percent only for the 
first and/or second modes, and these roughly parallel the errors in 
the Lin et. al . analyses. 
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Plots of chosen NASTRAN-computed eigenmodes are S res J n ™ 

Fiaures 3-8. Plate 761X does not appear since its mode shapes are 
identical to plate 761L, and plates 764, 770 and 772 are omitted 
since their mode shapes are very similar to those of plates 7j 4 r 
769 and 771. The nodal patterns can be confirmed as being identi- 
ca!'to?hose measured. (Figures 3-8 should be compared to Tables 8 
7 io 9 11 and 3, respectively, of reference 11). Finer details 

of some NASTRAN modes, particularly theveeringoinodallinesaway 
frnTn each other in plates with angled plies, are more easily seen m 
co?Srld grapMcs terminal displays, only the plates with specially 

orthotropio ply layups (734, 761L, 761X ' e *? ate sidas 

modes with nodal lines more or less parallel to the plate sides. 

Most of these modes are essentially beam-like fle^al ast^plates 
without phase changes at one symmetry plane. In contrast, pla 
with angled plies possess a larger number of more complex plate 

flexural modes. 

The natural frequency results obtained via NASTRAN CQUAD4 plate 
element idealizations in these simple composite plate vibration 
problems are judged to be acceptably accurate for engineering pur 
TheVlemint performs as well as alternative or.a a aons 
similar accuracy (the Lin et. al. element) for both GFRP and highly 
anisotropic CFRP plates for a variety of ply geometries. Potential 
usels ofihe element must be cautioned that this validation effort 
concerned plates with L/t ratios in the vicinity of 100 to *50. ^ or 
such L/t, modeling of laminate transverse shear stiffness helps to 
eke out 4 few percent in low-order natural frequency accuracy, but 
?s not absolutely essential to obtain rough-cut results. A more 
critical test of the CQUAD4 for modeling highly anisotropic lami 
nSei would have to concern plates with L/t lower than, say, about 

50. 

It should be mentioned that Lin et. al. and Xiao et. al. also 
measured and computed the specific damping capacities (2* times the 
modal loss factor) of each mode. Their FEA program was capable of 
modeling orthotropic lamina damping properties. No attempt “as made 
to predict modal damping factors in this effort, as NASTRAN is curr 
ently restricted to the modeling of isotropic material damping. 

SUMMARY 

The performance of the NASTRAN CQUAD4 membrane and plate 
element in analysis of undamped natural vibration modes of thin 
fiber-reinforced composite plates has been evaluated. The element^ 
provides natural frequency estimates that are comparable in accura 
cy to alternative formulations, and, in most cases, deviate by les 
than ten percent from experimentally measured frequencies. The 
predictions^ ie within roughly equal accuracy bounds for the two 
llltitll types treated (GFRP and CFRP) , and for the ply layups 
considered (unidirectional, cross-ply, angle-ply). Effective 

elastic lamina moduli had to be adjusted for m ® as " r ®^^ e at V f iftv 
fraction to attain this level of accuracy; nominal moduli at fifty 
percent volume fraction gave significantly inferior frequency 
estimates. The lumped mass option provided more accurate frequen- 
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cies than the consistent mass option. 


Th i s evaluation concerned only plates with L/t ratios on the 
order of 100 to 150. since the CQUAD4 utilizes first-order correc- 
tions for transverse laminate shear stiffness, the element should 
provide useful frequency estimates for plate-like structures with 
lower L/t. For plates with L/t below 20, consideration should be 
given to idealizing with 3-D solid elements. 


Based on the observation that natural frequencies and mode 
shapes are predicted with acceptable engineering accuracy, it is 
concluded that the CQUAD4 should be a useful and accurate element 
for transient shock and steady-state vibration analysis of Naval 
ship structures. 
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Table 1. Geometric and material parameters 
of ten fiber- reinforced plastic 
plates tested and analyzed by 
Lin et. al. and Xiao et. al. 


1 

Plate | Plate 
Ident. | Material 
Number | 

I 

| Plate 
| Length 
j L 

| (nm) 

Plate 

Thickness 

t 

(mm) 

| Fiber 
| Volune 
| Fraction 

' ..!* 

Mass 

Density 

<?3 
(kg/m ) 

Number 

of 

Plies 

Ply 

Layup 

| (degrees; see 
| Figure 2) 

Source 

of 

Data 

(ref. no.) 

734 

GFRP 

j 227.0 

2.05 

1 

1 0.451 

1813.9 

8 

| [0/90/0/901 

10,11,12 

761 L 

GFRP 

| 182.75 

1.64 

| 0.568 

1971.0 

8 

| [03 


761 X 

GFRP 

j 249.0 

2.28 

| 0.530 

1924.7 

8 

j [0] 

12 

765 

GFRP 

j 230.5 

1.45 

j 0.607 

2023.6 

8 

j [45/-45/45/-451 

1 ’I 

769 

GFRP 

j 224.2 

1.37 

| 0.621 

2041.7 

8 

j [0/90/45/-45] 

1 11 

771 

GFRP 

j 204.6 

2.11 

j 0.592 

2003.5 

12 

| [(0/-60/60) ] 

| 10,11 

762 

CFRP 

j 178.0 

1.58 

j 0.516 

1566.0 

8 

j [0] 

1 11 

764 

CFRP 

j 234.5 

2.12 

| 0.342 

1446.2 

8 

j [0/90/0/90] 

| 10.11 

770 

CFRP 

j 215.0 

1.62 

j 0.494 

1551.4 

8 

j [0/90/45/-45] 

| 11.12 

772 

CFRP 

| 215.6 

2.02 

j 0.618 

1636.4 

12 

j [<0/-60/60> 2 ] 

| 10.11.12 


NOTE: All laminates are symmetric about plate midplane. 


Table 2. Nominal in-plane effective elastic 
moduli of GFRP and CFRP laminae at 
fifty percent fiber volune fraction 


Material 


Fiber 

type 


Resin 

type 


“11 



l 

| (Gpa) 

1 

(ps i 

l 

(GPa) 

I 

(psj) 

l 

(GPa) 

l 

(psj) 1 



1 

l 

i 

1 

/10 ) 

l 


I 

/10 ) 

l 


l 

/10 ) | 


GFRP 

“Glass" | 

i 

DX210 | 37.78 

1 

5.48 

I 

10.90 

I 

1.58 

1 

l 

4.91 

l 

0.71 | 

0.3 


I 

epoxy | 

1 


l 


i 


I 


l 

1 


CFRP 

HM-S | 

DX210 j 172.7 

1 

25.0 

l 

7.20 

I 

1.04 

l 

3.76 

l 

0.55 | 

0.3 


1 

epoxy | 

1 


I 


i 


I 


I 

1 


DX210 


1 3.21 

1 

0.47 

I 

3.21 

I 

0.47 

l 

1.20 

I 

0.17 | 
| 

0.34 

epoxy 

1 1 

1 

1 


l 


I 


I 


I 

l 



22 


12 


V 


12 


NOTE: G = 0.6(G ) = 2.94 GPa is assuned for GFRP, 

G 23 = 0.2(E 12 ) = 1.44 GPa is assuned for CFRP. 
23 22 
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Table 3. In-plane and out-of-plane effective elastic 
moduli for GFRP and CFRP laminae, adjusted 
for fiber volume fraction of each plate 



Plate 

Plate 

Fiber 

e 

E , 



G 



I dent. 

Material 

Volune 


22' 

E 


G 12 ' 

23 



Number 


Fraction 

(GPa) 

33 


13 

(GPa) 





v f 


(GPa) 


(GPa) 




734 

GFRP 

0.451 



34.4 

9.7 


4.3 

2.6 



761 L 

GFRP 

0.568 

42.5 

13.0 


5.7 

3.4 



761 X 

GFRP 

0.530 

39.9 

12.0 


5.2 

3.1 



765 

GFRP 

0.607 

45.2 

14.2 


6.5 

3.9 



769 

j GFRP 

0.621 

46.2 

14.4 


6.7 

4.0 



771 

j GFRP 

0.592 

44.1 

14.0 


6.2 

3.7 



762 

| CFRP 

0.516 

178.0 

7.4 


3.9 

1.5 



764 

CFRP 

0.342 

119.0 

5.6 


2.6 

1.1 



770 

CFRP 

0.494 

171.0 

7.1 


3.7 

1.4 



772 

CFRP 

0.618 

213.0 

8.7 


5.0 

1.7 



NOTE: V 

^2 = 0,30 ,s assumed for all laminates, with 

no adjustment. 



Table 4. Comparison of natural frequencies from final 

NASTRAN CQUA04 idealizations and predictions of 
Lin and Xiao et. al. to measured data 


Plate 

Number 


Mode 

Number 


Natural Frequency Ratios, (computed / measured) 


and | 

pa rams . | 

i 


i 

1 

1 

l 

CLT 

ref. 12 

Lin et. | 
al. ref. | 
11,12 | 

Xiao 
et. al. 
ref. 12 

NASTRAN 

12-by-12 

lumped 

NASTRAN j 

12-by-12 j 
cons is. | 

734 j 

GFRP j 

1 

1 

l 

l 

1.20 

1.07 j 
| 

0.98 

1.03 

1.05 j 

8 plies | 

| 

2 

1 

I 

1.09 

1.00 j 

I 

1.00 

0.94 

0.97 | 

1 

10/90/0/ | 
90]sym | 

3 

1 

l 

1 

1.20 

1 

1.03 | 

1 

1.00 

0.98 

1.01 | 

1 

L/t = | 

4 

1 

1 

1.11 

1.05 | 

| 

1.03 

0.99 

1.03 j 

110.7 j 
1 

5 

1 

i 

1.12 

1.04 j 
| 

1.01 

1.00 j 

1.03 j 

i 

1 

| 

6 

1 

1 

1 

1.11 



1 

1.06 | 

1 

i- 

1.01 

0.98 | 

1 

1 

1.04 | 
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Table 4. (Continued) 


Plate | Mode | Natural Frequency Ratios, (computed / measured) 
Number I Number | 

an d I I CLT | Lin | Xiao | NASTRAN | NASTRAN 

l ' . I - . m ** I « ** ■ -in 


1 1 


I 

ref. 12 

1 

1 

ref. 11 

1 

ref. 12 

lumped 

consis. 

1 

1 1" 

I 761 L | 

1 

I-- 

I 

— 

i 

1 

1.13 


— 

1.11 

1.13 


| GFRP | 


I 


1 







| 8 plies | 

2 

l 

i 

— 

1 

i 

1.00 


— 

0.97 

1.00 


| [03 sym | 

3 

i 

I 

i 

— 

1 

1 

t 

1.05 


— 

1.02 

1.05 


i L/t = i 

4 

i 

I 

— 

i 

l 

1.00 


— 

0.94 

0.97 


| 111.4 1 


i 


l 








5 

I 

— 

1 

1.04 


— 

0.99 

1.02 



6 

I 

l 

I 

— 

1 

1 

1 

1.03 


— 

0.97 

1.03 


| 761 X | 

1 

I 

1.18 

1 

1.09 


0.98 

1.10 

1.11 


| GFRP | 


I 


1 







| 8 plies | 

2 

i 

1.20 

1 

1 

1.08 


1.00 

1.02 

1.05 


| [0] sym | 

3 

i 

I 

i 

1.20 

i 

1 

i 

1.08 


0.98 

1.05 

1.08 


| L/t = | 

4 

1 

I 

1.08 

i 

l 

0.99 


0.99 

0.94 

0.96 


j 109.2 | 


I 


1 








5 

I 

1.10 

1 

i 

1.02 


0.99 

0.98 

1.02 



6 

I 

i 

I 

1 

1.14 

1 

1 

1 

1.03 


0.99 

0.95 

1.01 


| 765 j 

1 

I 

— 

1 

1.09 


— 

0.82 

0.83 


| GFRP | 


I 


1 







j 8 plies | 

2 

I 

i 

— 

1 

| 

0.88 


- — 

0.65 

0.67 


| [4S/-45/ | 

3 

1 

I 

— 

1 

1 

0.95 


— 

0.68 

0.70 


| 45/-45] | 


I 


1 







| sym | 

4 

i 

i 

— 

l 

| 

1.07 


... - 

0.78 

0.81 


I L/t = | 

5 

1 

I 

— 

1 

l 

1.06 


— 

0.78 

0.81 


| 159.0 | 


I 


l 







l l 

l l 

i i 

6 

I 

i 

--i 


1 

1 

--i 

1.02 



0.72 

0.77 
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Table 4. (Continued) 


1 1 1 

| Plate | Mode | 

| Number | Number | 

| and j 

| pa rams. | | 

l 1 1 

Natural 

CLT | 

1 

ref. 12 | 

1 

Frequency Ratios, (computed / measured) 

Li " | Xiao | NASTRAN | NASTRAN 

et. al. j et. al. j 12-by-12 j 12-by-12 

ref. 11 | ref. 12 | lumped | consis. 

1 1 . 

769 | 

GFRP | 

1 

i 

1 

l 

— 

1 

1 

1 

0.95 

1 

r 

1 

| 

0.93 

" r 

0.95 

8 pi ies | 
1 

2 

1 

1 

— 

1 

1.08 

— 

1 

1.03 


1.06 

1 

[0/90/45/| 
-45] sym | 

3 

1 

1 

l 

— 

1 

1 

1 

0.98 

— 

1 

1 

1 

0.92 


0.94 

I 

L/t = | 

4 

l 

1 

— 

1 

1 

1.01 

— 

1 

| 

0.96 


0.99 

163.6 j 
1 

5 

1 

i 

— 

1 

1 

1.04 

— 

1 

l 

0.99 


1.03 

i 

l 

l 

I-. 

6 

l 

l 

— 

1 

1 

1 

1.01 

— 

1 

1 

1 

0.94 


0.99 

771 1 

GFRP | 

1 

1 

1 

1 

— 

1 

1 

1.20 

— 

1 

| 

1.19 


1.20 

12 plies | 

2 

l 

1 

— 

1 

1.17 

— 

1 

1.13 


1.16 

i 

[0/ - 60/60 1 
/0/- 60/60) | 

3 

1 

1 

1 

— 

1 

1 

1 

0.98 

— 

1 

1 

I 

0.93 


0.96 

sym | 

1 

4 

1 

i 

— 

1 

1 

1.06 

— 

1 

1 

1.02 


1.06 

L/t = j 

97.0 | 

5 

I 

1 

1 

— 

1 

1 

1 

1.07 

— 

1 

1 

| 

1.03 


1.07 

1 

1 



6 

l 

1 


1 

1 

1.03 

— 

1 

1 

_ _ 1 _ 

0.95 


1.01 

1 

762 | 

CFRP | 

1 

’ 1 ” 
1 
l 

— 

" 1" 
1 
1 

1.03 

— 

1 

| 

1.01 


1.02 

8 plies | 
1 

2 

l 

i 

— 

1 

1 

1.10 

— 

1 

l 

1.02 


1.05 

i 

[0] sym | 
1 

3 

1 

1 

1 

— 

1 

1 

I 

1.06 

— 

1 

1 

1 

0.99 


1.04 

L/t = | 

112.7 | 

4 

1 

1 

1 

— 

1 

1 

1 

1.11 

— 

1 

1 

| 

1.01 


1.08 

l 

i 

5 

l 

i 

— 

1 

1 

1.10 

— 

1 

1 

1.00 


1.07 

i 

l 

l 

I— 

6 

1 

l 

1 

1— 

— 

1 

1 

1 

-I-- 

1.03 

— 

1 

1 

1 

-i-- 

1.00 

i 

i 

i-- 

1.03 
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Table 4. (Continued) 


Plate | Mode | Natural Frequency Ratios, (confuted / measured) 
Member | Number | 
and | | 

pa rams. 


CLT 

ref. 12 

Lin 

et. al. 
ref. 11 

Xiao 
et. al. 
ref. 12 

| NASTRAN 
| 12-by-12 
| lumped 

NASTRAN 

12-by-12 

consis. 

— 

0.84 

— 

| 0.81 

| 0.82 

— 

0.97 

— 

j 0.95 

j 0.97 

— 

0.97 

— 

j 0.93 

| 0.96 

— 

0.99 

— 

| 0.97 

| 0.99 

— 

1.00 

— 

| 0.97 

j 1.00 

— 

0.98 



— 

j 0.92 

| 0.97 

1.14 

1.11 

1.08 

| 1.10 

j 1.12 

1.16 

1.11 

1.06 

| 1.09 

| 1.12 

1.15 

1.09 

1.00 

| 1.06 

j 1.10 

1.07 

1.00 

1.00 

| 0.98 

| 1.01 

1.15 

1.08 

0.99 

| 1.06 

| 1.09 

1.11 

1.03 

1.00 

j 0.99 

j 1.04 

1.07 

1.05 

1.00 

j 1.04 

j 1.05 

1.05 

1.03 

1.00 

j 0.99 

| 1.02 

1.07 

1.04 

1.01 

| 1.01 

| 1.04 

1.10 

1.06 

1.00 

j 1.02 

j 1.06 

1.09 

1.05 

1.00 

j 1.02 

j 1.05 

1.09 

1.03 

0.99 

| 0.97 

j 1.03 


764 
CFRP 
8 plies 

[0/90/0/ 
90) sym 

L/t * 

110.6 


1 

2 

3 

4 

5 

6 


770 
CFRP 
8 plies 

[0/90/45 

/-45]sym 

L/t = 
132.7 


1 

2 

3 

4 

5 

6 


772 

CFRP 

12 plies 

[0/ *60/60 
/0/-60/60] 
sym 

L/t = 
106.7 


1 

2 

3 

4 

5 

6 
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Table 5. Root -mean -square errors in NASTRAN - computed natural 
frequencies for ten composite plates for uncorrected 
and corrected effective lamina elastic moduli 


Plate 

Number 


RMS error in NASTRAN -predicted frequency 



1 

l 

With uncorrected 
moduli (percent) 

With corrected 
moduli (percent) 

734 

1 

l 

5.3 

2.9 

761 L 

1 

8.8 

5.5 

761 X 

1 

10.7 

5.7 

765 

l 

34.4 

26.7 

769 

l 

15.6 

5.4 

771 

1 

9.7 

10.1 

762 

l 

2.1 

0.8 

764 

l 

11.9 

9.6 

770 

1 

7.1 

6.6 

772 

l 

9.6 

2.5 


Table 6. Comparison of absolute measured and predicted 

composite plate natural frequencies of Lin et. al 
with NASTRAN C QUAD 4 computations 


Plate | Node 
Number I Number 


Natural Frequencies (Hz) 


and | 

pa rams. | 

l 


I 

l 

l 

Lin et. al., 
ref. 10-12 
measured | computed 

| NASTRAN 
| 12-by-12 
| lumped 

NASTRAN 

12-by-12 

consis. 

1 

734 | 

GFRP | 

1 

1 

I 

I 

1 

62.2 | 
I 

66.4 

| 64.2 

65.0 

6 plies | 
1 

2 

l 

| 

131.4 | 

1 

131.6 

| 123.8 

127.4 

1 

C0/90/0/ | 
90] sym | 

3 

1 

l 

I 

1 

159.2 | 

1 

164.5 

j 156.4 

160.9 

1 

L/t = | 

4 

l 

l 

180.5 | 

1 

189.8 

j 177.9 

184.0 

110.7 | 

1 

5 

l 

i 

200.1 | 
| 

208.9 

j 199.2 

206.0 

i 

l 

i 

6 

1 

l 

i 

1 

326.7 | 

1 

347.2 

1 321.2 | 

1 i 

338.4 

I 
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Table 6. (Continued) 


1 1 

| Plate | Mode 

| Number | Number 

I and j 

| pa rams. | 

I I 

■ i 

Natural Frequencies (Hz) 

Lin et. al., | NASTRAN | NASTRAN 

ref. 10-12 | 12-by-12 | 12-by-12 

measured | computed | lumped | consis. 

1 1 1 


1 I" 

| 761 L | 



1 

1 

78.1 | 

88.1 

86.9 

88.0 


| GFRP j 


1 





| 8 plies | 

2 

131.2 | 

130.7 

127.1 

130.7 


1 ! 

| [0] sym | 

■ i 

3 

211.5 

222.2 

215.0 

222.2 


1 ! 
| L/t = | 

4 

246.0 

246.1 

232.1 

238.7 


j 111.4 j 







l l 

5 

287.1 

297.8 

284.1 

294.2 


l 1 

l i 

l I 

6 

362.6 

374.4 

352.9 

375.0 


j 761 X j 

1 

57.2 

62.5 

62.9 

63.7 


| GFRP | 







| 8 plies | 

2 

90.3 

97.4 

92.5 

95.2 


1 1 
| [0] sym | 

3 

148.7 

160.5 

155.9 

161.2 


l l 

| L/t = | 

4 

181.6 

180.2 

170.4 

175.2 


j 109.2 | 







I I 

5 

211.2 

215.9 

207.5 

214.9 


I i 

I I 

I I 

• l 

6 

270.5 

278.8 

256.9 

272.9 


1 765 j 

1 

84.0 

91.3 

68.7 

69.6 


j GFRP j 







| 8 plies | 

2 

114.0 

99.9 

74.4 

76.6 


1 i 

| [45/ -45/ | 

3 

157.0 

149.5 

107.1 

110.0 


| 45/ - 45] | 







| sym | 

i I 

4 

199.3 

212.6 

156.3 

161.6 


i i 

| L/t = | 

5 

213.4 

226.7 

167.2 

173.4 


j 159.0 j 







1 1 

1 1 

* 1 

6 

346.6 

i 

353.5 

i 

249.6 

265.5 

--I 
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Table 6. (Continued) 


| Plate | 
| Number | 

l and i 

| params. | 

I | 

Mode | 

Number | 

Natural Frequencies (Hz) 

Lin et. al., | NASTRAN | NASTRAN 

ref. 10-12 j 12-by-12 | 1Z-by-12 

measured | computed | lumped | consis. 

1 1 1 . 

1 ! 
1 769 1 

1 

] 

1 

58.2 


1 

54.4 

I 

55.1 j 

| GFRP | 


l 


l 


i 


j 8 pi ies | 
1 1 

2 

l 

| 

91.6 

| 99.0 

1 

94.5 

i 

97.2 | 

1 1 
| [0/90/45/| 

3 

1 

1 

125.5 

1 

| 123.0 

115.2 

l 

1 

118.4 j 

| -45] sym | 


1 


l 


l 


l 1 

4 

1 

150.4 

| 151.3 

143.8 

l 

148.7 | 

| L/t = | 


l 


I 


i 


| 163.6 | 

| j 

5 

1 

1 

156.8 

| 163.0 

| 

155.7 

l 

161.1 j 

i 1 

l 1 

6 

1 

l 

l 

277.3 

I 

| 279.4 

1 

260.5 

i 

i 

I 

1 

274.4 | 

1 77 1 | 

1 

1 

90.4 

1 

| 108.2 

107.4 

l 

108.8 j 

1 GFRP i 


l 


1 


I 


| 12 plies | 
1 1 

2 

1 

j 

144.7 

| 168.6 
i 

163.6 

i 

168.4 j 

1 ! 
| [0/ - 60/60 | 

3 

1 

1 

222.3 

1 

| 218.6 

206.6 

I 

i 

212.4 j 

| /0/-60/60] | 


l 


i 


I 


I sym | 

| 1 

4 

l 

i 

264.1 

| 280.2 
| 

270.6 

I 

i 

279.8 j 

j L/t = | 

5 

1 

1 

281.1 

I 

| 301.0 

290.9 

i 

I 

300.9 | 

j 97.0 | 


1 


l 


l 


1 1 

1 l 

i I-. 

6 

1 

1 

492.6 

| 505.2 

| 

469.7 

l 

l 

1 

499.3 | 

1 76 ? | 

1 

" 1" 
l 

81.5 

| 83.6 

82.2 

I 

83.3 j 

| CFRP | 


l 


1 


l 


| 8 plies | 

1 [ 

2 

l 

i 

107.4 

| 118.4 

| 

109.4 

i 

112.5 | 

| [0] sym | 

| | 

3 

1 

1 

1 

196.6 

1 

| 207.8 

1 

198.4 

i 

l 

205.1 j 

1 L/t = | 

4 

1 

1 

295.5 

1 

| 329.4 

299.2 

I 

1 

318.2 | 

1 112-7 | 


1 


1 


l 


1 1 
1 1 

5 

1 

i 

382.5 

| 419.8 

| 

383.6 

I 

t 

410.2 | 

1 1 
1 1 

6 

1 

1 

531.0 

1 

| 546.9 

530.6 

i 

I 

545.9 | 
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Table 6. (Continued) 


Plate | Mode j Natural Frequencies (Hz) | 


Number | Number | 

and j I 

pa rams. | | 

i 1 

i 1 

Lin et. al. ( 
ref. 10-12 

measured | computed 

i 

| NASTRAN | NASTRAN 
j 12-by-12 | 12-by-12 
| Imped | consis. 

1 — 

764 | 

1 

1 

68.9 

i 

l 

58.1 

| 55.5 

| 56.2 

CFRP | 

8 plies | 

2 

218.9 


213.3 

1 

| 206.9 

I 

1 

| 212.8 
j 

1 

[0/90/0/ | 

3 

251.2 


243.5 

| 233.3 

| 241.5 

i 

90] syra | 

i 

4 

305.4 


302.5 

i 

| 294.9 

1 

| 303.4 

L/t = | 

110.6 | 

5 

323.5 


324.2 

1 

| 312.7 

1 

1 

| 323.7 

1 

1 

1 

1 

6 

452.5 


441.6 

1 

| 414.5 

1 

1 

| 437.8 

l 

I-- 

770 | 

1 

77.8 


86.3 

| 85.9 

1 

| 87.0 

CFRP | 

8 plies | 

2 

202.7 


224.5 

1 

| 220.6 
1 

1 

| 226.9 

| 

i 

[0/90/45 | 

3 

258.0 


280.4 

1 

| 273.9 

i 

| 283.4 

/-451 sym | 

I 

4 

298.7 


298.8 

l 

| 293.6 

1 

| 302.0 

L/t = | 

132.7 | 

5 

322.0 


348.4 

1 

| 340.6 

| 

1 

| 352.4 

j 

1 

1 

1 

6 

496.7 


512.2 

1 

| 490.0 

| 

| 517.3 

1 

772 j 

1 

156.6 


165.2 

1 

| 162.7 

l 

| 165.0 

CFRP | 

12 plies | 

i 

2 

272.0 


279.1 

1 

| 269.9 

1 

1 

| 277.6 

| 

i 

[0/ -60/60 1 

3 

372.3 


387.8 

1 

| 376.9 

| 387.6 

/0/- 60/60] j 

sym | 

4 

407.8 


432.6 

1 

| 416.7 

| 

1 

| 431.1 

| 

1 

L/t = | 

5 

486.1 


511.4 

1 

| 494.7 

| 512.0 

106.7 | 

1 

1- 

6 

779.0 

I 

_l 

800.4 

1 

| 752.6 

1 

| 799.4 

1 


I 

I 

I 
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Table 7. Mode-by-mode tabulation of NASTRAN CQUAD4 idealization 
frequency prediction errors, with averaged percentage 
errors and eigenmode descriptions 


Plate 

Mode 

i 

( Percentage errors 

1 

| Eigenmode 

Number 

Number | 



j description 

and 


| Lin et. 

| NASTRAN 

| NASTRAN 

i 

params. 


1 al. 

| 12-by-12 | 12-by-12 1 

I 

1 


1 

| lumped 
.. 1 

| consis. 

i 

i 

i 

l 

734 | 

1 

1 

1 ♦ 7 

i 

1 ♦ 3 

i 

1 ♦ 5 

| shear and flexure 

GFRP | 


1 

1 

| 

i 

8 plies | 

2 

1 o 

1 ' 6 

1 ' 3 

| 2-noded beam abt. 



1 

1 

1 

| weak ax i s 

[0/90/0/ | 

3 

1 ♦ 3 

1 " 2 

1 + 1 

| 2-noded beam abt. 

90] sym | 


1 

i 


| strong axis 


4 

1 ♦ 5 

1 ' 1 

1 + 3 

| 2-noded beam abt. 

L/t = j 


I 

I 


| weak axis + shear 

110.7 j 

5 

1 + 4 

1 0 

1 + 3 

| 2-noded beam abt. 



i 

I 


| strong axis + shear 


6 

1 + 6 

1 " 2 

1 + 4 

| fundamental 



i 

j 


| plate flexure 


avg. 

| 

| + 4.2 

1 " 1-3 

j + 2.2 


761 L j 

1 

1 ♦ 13 

! ♦ ii 

1 ♦ 13 

| shear and flexure 

GFRP | 


i 

i 



8 plies | 

2 

1 0 

1 • 3 

1 o 

| 2-noded beam abt. 



I 

l 


| weak axis 

[0] sym j 

3 

1 + 5 

1 + 2 

1 + 5 

| 2-noded beam abt. 



i 

i 


| weak axis ♦ shear 

L/t = j 

4 

1 o 

1 ’ 6 

1 * 3 

| 2-noded beam abt. 

111.4 | 


l 

1 





1 

1 


| strong axis 


5 

1 + 4 

1 * 1 

1 + 2 

| 2-noded beam abt. 



1 

l 


| strong axis + shear 


6 

1 + 3 

1 - 3 

1 ♦ 3 

| 3-noded beam abt. 



1 

1 


f weak ax i s 


avg. 

| + 4.2 

1 0 

j + 3.3 


761 X | 

1 

! + 9 

j + 10 

1 + 11 

| shear and flexure 

GFRP | 


1 

1 



8 plies | 

2 

1 ♦ 8 

1 + 2 

1 ♦ 5 

| 2-noded beam abt. 

1 


1 

1 


| weak axis 

CO] sym | 

3 

1 + 8 

1 ♦ 5 

j + 8 

| 2-noded beam abt. 

1 


1 

1 


| weak axis + shear 

L/t = | 

4 

1 - 1 

1 ' 8 

j * 4 

| 2-noded beam abt. 

109.2 | 


1 

1 


| strong axis 

1 

I 

5 

1 + 2 

1 * 2 

1 + 2 

| 2-noded beam abt. 

1 


1 

1 

1 

| strong axis + shear 

1 

6 

1 + 3 

1 - 5 

1 > 1 

1 3-noded beam abt. 

1 

I-- 


1 

1 

1 

1 

1 

i 

| weak axis 

1 

1” 

avg. 

1 

| ♦ 4.8 

1 

1 

| + 0.7 

1 

i 

| + 3.8 
1 

i i 

i i 


76 


Table 7. (Continued) 


Plate | 
Number | 
and | 

pa rams. | 

Mode 

Number 

Percentage errors 

Lin et. | NASTRAN | NASTRAN 

al. | 12-by-12 j 12-by-12 

| lumped | consis. 

i i 

| E i genmode 

| description 

765 j 

1 

+ 9 

| - 18 

1 

1 • 17 

| plate flexural 

GFRP j 



l 



6 plies | 

2 

- 12 

1 ‘ 35 
i 

1 - 33 

| plate flexural 

[45/ -45/ | 

3 

- 5 

1 

1 “ 32 

| - 30 

| plate flexural 

45/-45] j 



I 



sym | 

4 

♦ 7 

| - 22 

1 

j - 19 

| plate flexural 

L/t = j 

5 

♦ 6 

1 

1 " 22 

j - 19 

| plate flexural 

159.0 j 



l 




6 

♦ 2 

| - 28 
1 

1 

1 ' 23 

| plate flexural 


avg. 

+ 1.2 

1 

| - 26.2 
1 

j - 23.5 


769 j 

1 

- 5 

1 

1 ' 7 

1 " 5 

| shear and flexure 

GFRP | 



1 



8 plies | 

2 

+ 8 

1 + 3 

1 + 6 

| 2-noded beam abt. 




l 


| weak ax i s 

[0/90/45/ j 

3 

- 2 

1 “ 8 

1 ' 6 

| 2-noded beam abt. 

-45] sym | 



1 


| strong axis 


4 

♦ 1 

1 ' 4 

j - 1 

| plate flexural 

L/t = j 



l 



163.6 j 

5 

+ 4 

1 * 1 

1 + 3 

| plate flexural 


6 

♦ 1 

1 

1 “ 6 
1 

i _ 

1 " 1 

j plate flexural 


avg. 

♦ 1.2 

I 

| - 3.8 

1 

| - 0.7 


771 | 

1 

+ 20 

1 

| + 19 

j + 20 

| shear and flexure 

GFRP | 



1 



12 plies | 

2 

+ 17 

1 + 13 

j + 16 

| 2-noded beam abt. 




i 


| weak axis 

[0/-60/60] 

3 

- 2 

1 * 7 

1 * 4 

| 2-noded beam abt. 

/0/-60/60] | 



1 


| strong axis 

sym | 

4 

+ 6 

1 + 2 
1 

1 + 6 

| plate flexural 

L/t = j 

5 

+ 7 

1 

1 + 3 

1 ♦ 7 

| plate flexural 

97.0 j 



1 




6 

♦ 3 

1 " 5 

1 + 1 

| 3-noded beam abt. 




1 

1 


| weak ax i s 

1 

avg. 

♦ 8.5 

1 

1 

| ♦ 4.2 

1 

1 * 7 - 7 
1 

1 
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Table 7. (Continued) 


Plate | 

Mode 

i 

| Percentage errors 

| Eigenmode 

Number ) 

Number 




[ description 

and | 


| Lin et. 

| NASTRAN 

| NASTRAN 


pa rams. | 


| al. 

j 12- by- 12 

j 12-by-12 





| lumped 

| consis. 


762 j 

1 

1 + 3 

1 4 1 

| ♦ 2 

| shear and flexure 

CFRP | 



l 



8 pi ies | 

2 

j + 10 

1 4 2 

1 4 5 

j 2-noded beam abt. 




i 


| weak axis 

[0] sym | 

3 

1 4 6 

1 " 1 

j + 4 

| 2-noded beam abt. 




I 


| weak axis + shear 

L/t = j 

4 

1 4 11 

1 4 1 

j ♦ 8 

| 3-noded beam abt. 

112.7 j 



i 


| weak axis 


5 

j ♦ 10 

1 0 

j + 7 

| 3-noded beam abt. 




i 


| weak axis + shear 


6 

j + 3 

1 0 

1 4 3 

| 2-noded beam abt. 




i 


| strong axis 


avg. 

1 4 7 * 2 

j + 0.5 

| + 4.8 


764 j 

1 

j - 16 

1 ’ 19 

| - 18 

| shear and flexure 

CFRP | 



I 



8 pi ies | 

2 

1 * 3 

1 ‘ 5 

| - 3 

| 2-noded beam abt. 




i 


| weak ax i s 

[0/90/0/ j 

3 

1 ' 3 

1 - 7 

j - 4 

| 2-noded beam abt. 

90) sym | 



l 


| weak axis + shear 


4 

1 " 1 

l • 3 

1 • 1 

| 2-noded beam abt. 

L/t = j 



I 


| strong axis 

110.6 j 

5 

1 0 

l • 3 

1 0 

| 2-noded beam abt. 




l 


| strong axis + shear 


6 

1 " 2 

l • 8 

j - 3 

| fundamental 




l 

1 


| plate flexure 


avg. 

j - 4.2 

1 

1 " 7 - 5 
1 

| - 4.8 


770 j 

1 

j + 11 

1 

| ♦ 10 

1 4 12 

| shear and flexure 

CFRP | 



1 



8 pi ies | 

2 

1 4 11 

1 4 9 

1 4 12 

| 2-noded beam abt. 




1 


{ weak axis 

[0/90/45 j 

3 

1 4 9 

1 4 6 

| ♦ 10 

| plate flexural 

/-45]sym | 



1 




4 

1 0 

1 - 2 

1 4 1 

| 2-noded beam abt. 

L/t = j 



1 


| strong axis 

132.7 j 

5 

j ♦ 8 

1 4 6 
j 

| + 9 

| plate flexural 


6 

1 4 3 

1 

1 ' 1 
1 

1 

1 4 4 

| plate flexural 

1 

avg. 

j + 7.0 
1 

1 

| + 4.7 

1 

j + 8.0 
1 

1 
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Table 7. (Continued) 


Plate 

Number 

and 

pa rams. 

Mode 

Number 

Percentage errors 

Lin et. | NASTRAN | NASTRAN 

al. | 12-by-12 | 12-by-12 

| lumped | consis. 

i i 

| E i genmode 

| description 

i 

i 

| 

772 

1 

+ 5 

1 I 

1 + 4 

♦ 5 

I 

| shear and flexure 

CFRP 



I 


i 

12 plies 

2 

♦ 3 

1 ' 1 

♦ 2 

| 2-noded beam abt. 




l 


| weak axis 

[0/-60/60 

3 

+ 4 

1 + 1 

♦ 4 

| 2-noded beam abt. 

/0/-60/60] 



I 


| strong axis 

sym 

4 

♦ 6 

1 * 2 

♦ 6 

| plate flexural 
| 

L/t = 

5 

+ 5 

1 

1 * 2 

♦ 5 

1 

| plate flexural 

106.7 



I 


i 


6 

♦ 3 

I • 3 

+ 3 

| 3-noded beam abt. 




l 


| weak axis 


avg. 

+ 4.3 

| + 0.8 

+ 4.2 

' 1 
i 
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Notes : Positive X-axis 
Positive Y-axis 


fiber orientation angle of 0 degrees 
fiber orientation angle of 90 degrees 


8-noded rectangles 4-noded CQUAD4's 




(i) Lin and Xiao et. al. 
idealization 


(ii) NASTRAN CQUAD4 
idealization 


Ficrure 2. Plate element meshes used by Lin and 
Xiao et. al. and in NASTRAN analyses 
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PUflTTE tO. VIA OF XIRO.LIN.JU C 0/90/0/90 3S CFRP < |g -O V - I8> 
SOLUTION I - BIGB+JPUJE WhRLVSIS 
CSC c:r*4 9 .EIGU I .065 1 E+OC3 


Fig. 3c. Plate 734, Mode 3 





Fig. 3f . Plate 734, Mode 6 

Figure 3. (Continued) 
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PLflTE T» OF 
SOLUTION 3 - 

esc 7;hiN 


UIN,HI,RDW« <I9B*> 

eigemuwjle «r«-vsis 
T;EIQY a.S6^3E+GS3 


C C 0)03 GFWP 


Oa-BV-l2> 


Fig. 4a. 


Plate 761L, Mode 1 



Fig 4b. Plate 761L, Mode 2 


Figure 4 . 


NASTRAN-computed eigenmodes for plate 761L 




SOLUTION S - EIGENUffLUE flWLVSIB 
CSC li;W ll.'EIGU 8.79606+0fi3 



Fig. 4e. 


Plate 761L, Mode 5 



Fig. 4f. Plate 761L, Mode 6 

Figure 4. (Continued) 
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Figure 5. (Continued) 
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Fig. 5f. Plate 765, Mode 6 

Figure 5. (Continued) 
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PLHTE 709 OF Lffa.NI .RCFTS C I904> C0.90.4S, -4S3S GFRP < le-^-ia) 
SOLUTION 3 - EIGENUHLUE RNHLV5IS 
CSC 7;W 7;EIQU B . 77426+043 


Fig. 6a. Plate 769, Mode 1 



SOLUTION a - EIGBNURLUE ATRLSS1S 

esc B;m 6;EioJ e.BOBTE+osa 

Fig 6b. Plate 769, Mode 2 

Figure 6. NASTRAN-computed eigenmodes for plate 769 


9 
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CSC e;m Q.EICV e.7SME+0S3 

Fig 7b. p; 

Figure 7 . NASTRAN-compi 


9 



JCH 

CSC 


VIB PLATE TT 1 OP LIN#NI # FORMS C CO, -€0. C0>S3S GFHP < tg-gi*-te> 
SOUiTlON ] - EIGBMUHUUE FTWLV5IS 
lO.'MN 10 ;EICW e.3C7SE+OC3 


Fig. 7d. Plate 771, Mode 4 

Figure 7. (Continued) 
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Fig. 7e. Plate 771, Mode 5 



esc ib:w ta.Eicy 7. i4*ie+o« 

Fig. If. Plate 771, Mode 6 

Figure 7. (Continued) 
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Fig. 8a. Plate 762, Mode 1 



Fig 8b. Plate 762, Mode 2 


Figure 8 . 


NASTRAN-computed eigenmodes for plate 762 
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JCM VI 8 PUTTE 765 OF LIN.NI.nORM5 C(0>83 CF WP (15-BV-12) 
SOLUTION 3 - E I GENVRLUE RNRLVSIS 


Fig. 8c. Plate 762, Mode 3 



Fig. 8d. Plate 762, Mode 4 

Figure 8. (Continued) 
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Fig. 8e. Plate 762, Mode 5 



Fig. 8f. 


Plate 762, Mode 6 


Figure 8. 


(Continued) 
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